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We introduce a variant to the Diffusion Monte Carlo algorithm that can be employed to study the 
effects of the Rashba interaction in a many electron systems. Because of the spin-orbit nature of 
Rashba interaction a standard algorithm cannot be applied and therefore a specific imaginary time 
spin dependent propagator has been developed and implemented following previous work developed 
in the framework of nuclear physics. We computed the ground state energy of the 2D electron gas 
at different densities for several values of the Rashba interaction strength as a function of "Rashba 
spin states" polarization. Comparison is given with analytically known Hartree-Fock results and for 
the system in absence of Coulomb interaction. 
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I. INTRODUCTION 



The Rashba interaction, an electric field induced spin- 
orbit coupling, has been experimentally observed in semi- 
conductor hetero-structures, depending on their sym- 
metry, and has been proved to be tunable in strength 
through a gate voltageii^. The external voltage can then 
be used to control the spin state of the system. This 
property becomes extremely interesting in view of spin- 
tronics applications. Several experiments have also been 
performed with the aim of studying the dependence of 
the interaction strength on applied gate voltage and well 
thickness^. 

The typical experimental setup consists of a device 
etched on a two-dimensional quantum well formed at the 
interface of two semiconductors. In such devices elec- 
trons are constrained to move in a two-dimensional space 
(Ox—Oy) forming a 2D gas. The asymmetry of the quan- 
tum well generates an electric field, along the z-direction, 
perpendicular to the plane containing the electrons. This 
causes electrons to be subject to an effective magnetic 
field Be// oc p X E coupling to their spin. Such coupling 
gives rise to the well known Rashba potential^: 
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where is the momentum of the i — th electron, and erf 
and af are the Pauli matrices acting over the spin of par- 
ticle i. Neglecting the Coulomb interaction among the 
electrons the Hamiltonian is a sum of one body terms, 
and the problem is analytically solvable. Single parti- 
cle solutions are given by plane waves with k-dependent 
spinors. In particular, for each wave vector with momcn- 
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These solutions correspond to two different spin states, 
with the following dispersion law: 



£(fc) = |l±A|k| 



(3) 



This spin splitting causes non interacting electrons to 
arrange in two bands with different fillings depending 
on the interaction strength A, therefore inducing a nat- 
ural polarization of the system. In this case two dif- 
ferent Fermi surfaces are generated, as also shown ex- 
perimentally by beating patterns in the Shubnikov-De 
Haas oscillationsii^. Being A proportional to the inten- 
sity of the external electric field, it is possible to manip- 
ulate the occupation numbers and to induce at the same 
time precession on the electron spin. These properties 
induced new studies in the direction of spin field effect 
transistors^. 

Once the Coulomb interaction is included in the Hamil- 
tonian, no analytic solution is available. Several methods 
can be used to study this system and much theoretical 
work has been done so far (for a review see e.gi^). The 
Hartree-Fock method gives a simple analytical solution, 
and though it totally ignores the effect of correlations, 
it provides a useful insight on the structure of the single 
particle levels. Another interesting approach consists in 
applying a unitary transformation giving, to leading 
order in the spin-orbit strength, a transformed Hamilto- 
nian H = U~^HU whose eigenstates are also spin and 
angular momentum eigenstates. Though approximate, 
this method allows the use of standard Quantum Monte 
Carlo techniques^. 

In this work a particular implementation of the 
Diffusion Monte Carlo, directly dealing with spin-orbit 
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interactions, is proposed in order to obtain an ab-initio 
method giving accurate ground state energies of the sys- 
tem. No transformation is needed because the spin-orbit 
term is included in the imaginary time propagator. 

The paper is organized as follows: in Sec. [TI] we de- 
scribe the Diffusion Monte Carlo method and the way 
in which the Green's function used for treating Rashba 
interaction is derived. The trial wave function used as a 
starting point for the projection is described in Sec. Ill CI 
Finally in Sec. Illll thc results arc presented and discussed 
together with Hartree-Fock energies for the ground state 
of the system, and conclusions are given in Sec. IIVI 



II. METHOD 

A. Diffusion Monte Carlo method for general 
spin— independent many— body problems 

The Diffusion Monte Carlo (DMC) method is based 
on imaginary time evolution projection. Through the 
use of an appropriate Green's function, an initial state 
is projected over its lowest energy component having the 
same nodes or the same phase as the trial wave function 
ipT- The propagation in imaginary time r is given by 

V'(R,t) = exph(if-i;o)T]V(R,0) = 

cn exp - Eo)t]M^) , (4) 

n 

where Eq is a normalization factor, and R represents 
the space coordinates of the system. The amplitudes of 
higher energy states, due to propagation, decay expo- 
nentially with T, with a lifetime inversely proportional to 
their energy relative to Eq. If '0(R, 0) is not orthogonal 
to the ground state, '0(R, t —^ oo) will be proportional 
to the ground state itself. 

A typical many-body system is described with a 
Hamiltonian of the form 



Go is the exact Green's function of a two-dimensional 
non-interacting system: 



H = T + V, 



(5) 



where T is the sum of single-particle kinetic energy op- 
erators, and the potential V between all the electrons 
only depends on spatial coordinates. The evolution in 
imaginary time can be achieved by means of the Green's 
Function of the Hamiltonian that can be approximated 
by using Trotter's formula: 

e-^^ = e-^^e-^^ + 0{t). (6) 

In space representation the propagator ^ can be written 
as: 



G(R,R',t) = (R|e-(^-^«)"|R') 
~e[^W-^°l^Go(R,R',r). 
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where D = /m is the diffusion constant. The exponen- 
tial of the kinetic term gives rise to a free particle imag- 
inary time Green's function Gq, while the other term is 
viewed as a weighting term. 

The algorithm makes use of walkers, i.e. points in the 
coordinate space, in order to sample the ground state 
wave function. The propagation is obtained by diffus- 
ing the walkers according the displacements distribution 
given by Gq. Afterwards a weight is assigned to the walk- 
ers according to the other factor in the propagator given 
by the potential energy of the system. Because of the 
Trotter approximation, the Green's function is correct 
only at order 0{t). The problem is overcome by using 
short time propagation repeatedly in order to achieve 
long enough imaginary times and adequate statistics, 
minimizing at the same time the time step error. The 
result is given after an extrapolation to At — > 0. 

In order to enhance the efficiency, importance sampling 
is implemented in DMC algorithm a^°i^^ . The idea con- 
sists in sampling a density of points proportional to the 
ground state distribution multiplied by an importance 
function. In standard cases the importance function will 
only depend on space coordinates being the spin fixed, 
thus: 

^T(R)</'(R,r) = 

This can be shown to introduce a drift in Gq and to 
modify the weighting factor, which will not only con- 
tain the potential, but also the local energy El = 
i7'0T(R)/'0T(R) of the system. 



B. Many— body spin— dependent Hamiltonian 

The system we are studying is a two-dimensional elec- 
tron gas at T = 0, in presence of both Coulomb and 
Rashba interaction, with the addition of a uniform charge 
backgroundi2ii^. The Hamiltonian for the system can be 
written as 



H 



T.^ + >^T.(p>^-p^''') + ^^ouiiR), (10) 

i=l i=l 



(7) 



where Vcoui includes the electron-electron interaction 
and the effects of the background. The Hamiltonian in 
this case does not only contain space coordinate depen- 
dent potentials. The Rashba interaction, because of its 
spin-orbit character, contains spin and momentum op- 
erators, which therefore cannot be treated like simple 
weighting factors as in the case shown above. In order 
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to apply the DMC technique, a new propagator form is 
needed, taking into account the particular features of non 
locality and spin-dependenc o^'^'^^ . 

In order to simplify the following calculations let us 
first consider the single-particle Hamiltonian given by the 
kinetic term and the Rashba interaction for only one elec- 
tron. By applying Trotter's formula the following form 
of the Green's function can be obtained 

G(r, r', At) = e-^(P«'"--P-'"«)'^^Go(r, r', At) , (11) 

where the Pauli matrices and ay act on the spin com- 
ponents of the electron, and r and r' are the coordinates 
of the electron after and before the diffusion generated 
by Gq. The last factor Go will again give the space dis- 
placement, while the second one, including the Rashba 
interaction, contains both momentum and spin opera- 
tors, and can be viewed as acting on the free propagator. 
Expanding to first order in At the first factor of Eq. (fTTj) . 
and applying the derivatives given by p = — i/iV to Gq, 
we have 



1 - i^ia^Ay - ay Ax) 



exp 



-i — {axAy - ay Ax) 



Go(r,r',AT)~ 
Go(r,r',AT). (12) 



The first part can be interpreted as a spin rotation de- 
pending on Ax ~ x~x' and Ay = y—y', i.e. the displace- 
ment generated by the Gaussian free particle Green's 
function Gq. The appearance of the spin-rotating factor 
in the Green's function implies that the spin coordinates 
of the electrons must be explicitly used. Making use of 
a spinorial representation, the spin-state for the i ^ th 
electron is given by 



(13) 



where a and /3 are the amplitudes of the spin state in the 
{It)} and {U)} basis. 

The propagation of the spin-dependent Green's func- 
tion can then be realized with a rotation of spinors. For 
each electron i the spin-dependent propagator can be 
written using the following matrix form: 



exp 



-i^{axAy - ayAx) 



cosf-^Ar,) sin(AAr 

Ari 



-iAyi+Axj 



X 



An 



sm [^Ar 



cos {^An 



(14) 



where Ar, = Axf + Ayf . Therefore, during the prop- 
agation in imaginary time, an electron with initial coor- 
dinates (ri, ttj, /?i), will be first moved to (r[,ai,l3i) due 
to the free propagator Go, and then to (r-, a^, /3^) due to 
rotation according to 



(15) 



The approximations introduced after Eq. (jlip gener- 
ate an error which must be now taken into account at 
least to order At. If the spin-orbit propagator just de- 
rived were correct, we would expect to obtain 



^-t^{a^Ay-ayAx)Q^(^^^ r', AT)^A(r')dr' = 
^-x(pya^-p^ay) /"Go(r,r',AT)?A(r')dr'. (16) 



This equation does not hold unless the propagator is cor- 
rected for a weighting factor 
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It is important to note that the results obtained fol- 
lowing this approach coincide with the analytical form 
of the Green's function for the system in absence of 
Coulomb potential. In this case no error is introduced in 
G(R, R', At) because the Rashba interaction commutes 
with the kinetic term of the Hamiltonian. 

By considering the full many-particle Hamiltonian, in- 
cluding the Coulomb interaction, the total Green's func- 
tion is 

G(R,R', At) = ,-[v'c„.a(R)-i^„-^]A.^ 

e-*^E£,(Arr<-A.r.r)A-Go(R,R', At) . (18) 

When a spin-orbit term is introduced in the Hamil- 
tonian, as shown before, the propagator will not leave 
spin unchanged because space coordinate diffusion be- 
comes related to spin rotation. This means that also the 
importance sampling should not be naively applied, but 
in the total weight, a factor '!/'t(R, <5')/V't(R', <5'')' where 
S' and S respectively are the old and new spin states, 
should be taken into account. It is possible to apply the 
importance sampling adding a drift term in Gq as usual 
by considering 



^^(R,^) _ 
V't(R',^') 
V>t(R ,^0 ^t(R,^) 



V't(R,5') Vt(R',^) 



Vt(R',^') V't(R,^') ^ V^t(R',5') V^t(R',^') 



(,19) 



where the two forms of Eq. (1191) ^^'C equivalent to first 
order in At. The factor where spins are unchanged can 
be included as usual, while the second one can be in- 
terpreted as an additional weighting factor. We point 
out that in this case the local energy appearing in the 
weight includes only the the spin-independent part of 
the Hamiltonian. 



C. Wavefunction 

We used a trial wave function made of a Slater deter- 
minant of single particle states multiplied by a Jastrow 
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factor accounting for correlations, according to the form 



III. RESULTS 



N 



^t(R,^)-^(R,5)cxp -^ii(lri-rjl) 



(20) 



where u{r) is a two body pseudopotential of the double 
Yukawa fornJ^ii^. 

Single particle states used to build the determinant are 
chosen as the eigenstates of the Hamiltonian in absence of 
Coulomb interaction, coinciding also with Hartree-Fock 
single particle solutions, written in Eq. ^ 



t>k,±(r. 



a, 13) 



±a- 



ikr 



-f3 



(21) 



where a and /3 are respectively the up and down spin 
components with respect to the z axis. These are pairs 
of wave functions which we will call quasi-up (plus sign) 
and quasi-down (minus sign) states. 

DMC simulations of infinite systems are usually per- 
formed in close shell configurations in order to obtain 
a real wave function and to reduce the finite size ef- 
fects mainly due to the kinetic energy. This makes 
it possible to apply the fixed node approximation in 
order to obtain a state with the same symmetry as 
tpT- In our case the trial wave function cannot be re- 
duced to a real form because of the phase change in- 
duced by spin rotation. Therefore we employ the fixed- 
phase approximatio n^^i^^i^^i^^i^^i^° . For any complex 
wave function it is always possible to factorizc the mod- 
ulus in the following way 

M'f<-,S) = \M^,S)\e''-^^''^ . (22) 

Using fixed node approximation means finding the low- 
est energy state 0(R) whose product with the (in 
this case real) trial wave function is positively defined 
0(R)V't(R) = 0(R)'0t(R) > 0. In our case this product 
does not have a defined sign and cannot be used as a 
sampling density: 

<^(R,5)Vt(R,^) = 

\cj){R, 5)||Vt(R, 5)1 exp[i(0(R, S) - 0t(R, 5))](23) 

In the fixed phase approximation the problem is over- 
come assuming 



exp[i(6'(R, S) - 9t{'R, S))] = 1 , 



(24) 



which corresponds to finding a solution with the same 
phase as the trial wave function. Under this assump- 
tion Eq. becomes positive definite and the sampled 
distribution will only depend on the wave functions mod- 
ulus. 

In order to calculate the ground state energy at differ- 
ent polarizations and to reduce finite size effects. Twist 
Averaged Boundary Conditions (TABC) have been in- 
troduced in the algorith m^^i^^ . As shown in the electron 
gas without spin-orbit interactions, this allows to obtain 
good results without the need to use a very large num- 
ber of electrons, giving therefore a great improvement in 
terms of the computational time required. 



In order to verify our method a check has been per- 
formed by comparing numerical results with the exact 
solution obtained neglecting the electron-electron inter- 
action. As mentioned in the introduction, the Hamilto- 
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can be diagonalized by the single particle wave functions 
PT|) . The analytical solution to N fermions only subject 
to Rasliba interaction is therefore a N x N Slater deter- 
minant of these single particle solutions. For a given A 
value, the ground state has a well defined polarization. 
As a test for our DMC we multiplied the exact wave func- 
tion by a Slater determinant in order to obtain a wrong 
state with the same nodes as the real ground state. Our 
algorithm indeed proved to be able to project out the 
ground state obtaining the analytical result within error- 
bars. 

In this work we make use of the dimensionless param- 
eter Ts = ro/ao, defined as a function of the Bohr radius 
flo = f? jme^ and the radius containing only one particle 
on average V/N ~ Trrp. Energies are given in units of 
Rydberg/electron (1 Ry = e^/2ao) and the Hamiltonian 
of Eq. pU]) in these units is 
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(26) 
(27) 



where Vback contains the effects of the charge background. 

The main numerical results presented in this work con- 
cern the ground state energies of the two-dimensional 
electron gas with Coulomb and Rashba interaction at 
different densities (expressed in terms of the Vs param- 
eter), different values of the intensity A of the Rashba 
interaction and different polarizations defined as 



N 



(28) 



where N- and iV+ are the numbers of electrons in the 
quasi-down and quasi-up spin states respectively, and 
the total number of electrons in the simulation is fixed 
to = 58. All the values obtained for the ground state 
energies as a function of r^, A and ^ are reported in the 
appendix, where the DMC and the Hartree-Fock results 
are compared. 

In particular the Ts values 1, 5, 10 and 20 were cho- 
sen and, for each density, calculations were performed 
at different A values of ground state energy as a func- 
tion of quasi-up - quasi-down occupation numbers, the 
total number of electrons in our simulation being fixed 
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to N = 58. Coulomb interaction was treated by means 
of the Ewald sums2^. As a comparison, besides DMC re- 
sults, we also give Hartree-Fock energies, which are given, 




Polarization 




Polarization 



FIG. 1: (color online) Energy values at = 1 A = 0.1 re- 
spectively obtained with DMC (triangles), HF (squares) and 
without Coulomb interaction (circles). 

DMC results at non zero polarization are obtained by 
projecting out of a trial wave function ipx whose Slater 
determinant contains different numbers of quasi — up and 
quasi — down spin states. The approach is apparently 
similar to what is commonly done for the electron gas in 
absence of spin-orbit interaction. However, in this case 
the Hamiltonian does not commute with the z spin com- 
ponent. The choice of the quasi-up, quasi-down basis 
was made by comparison with the results obtained us- 
ing the standard spin-up and spin-down basis. The lat- 
ter consistently gives higher values of the energy for the 
ground state polarization. Results in absence of Coulomb 
interaction predict spin state polarization as a conse- 
quence of a two band dispersion law. A very similar 
behavior for the energy as a function of polarization is 
found in our results (see Fig. [1] ) . This suggests the exis- 
tence of an analogous two band structure also in presence 
of Coulomb interaction as experimentally proved by the 
Shubnikov-De Haas oscillations. 

In Fig. [2] energies at different A are shown as a func- 
tion of the polarization ^. As expected, the polarization 



FIG. 2; (color online) Energy values obtained at rs = 1 
for three different A values, respectively 0.1 (triangles), 0.2 
(squares) and 0.5 (circles). The arrows indicate the point 
of minimum energy for each A. For increasing A the energy 
minimum shifts to larger (in modulus) polarization. 



depends on the intensity of the spin-orbit interaction, in- 
creasing in module with larger A. Calculations have been 
performed at different rs values, corresponding to differ- 
ent densities. The onset of polarization occurs at lower 
and lower strengths of the Rashba interaction when the 
density is decreased. This is an effect of the Rashba in- 
teraction dependence on the momentum. In the Hamil- 
tonian of Eq. (P7)) the l/rg prefactors for the interaction 
terms, compared to the l/r^ of the kinetic energy let 
interactions become more efficient at high (low den- 
sity) compared to the kinetic term. While at high density 
the system will look less interacting, and therefore also 
less polarized, at larger values of rs the system will be 
more sensitive to the the spin-orbit interaction and will 
be more polarized also for smaller A. 

At lower densities correlations among electrons play 
a major role. In Fig. [T] and [3] DMC energies are shown 
together with the corresponding results from HF approxi- 
mation. We see a worsening agreement when increases. 
This must be due to the fact that HF approximation does 
not take correlations into account, causing a larger de- 
viation when the N-electrons wave function becomes less 
similar to a Slater determinant. 

In Fig. m energies for the 2D electron gas from Refi^ 
are compared with ground state results in presence of 
spin-orbit interaction with fixed strength A = 0.1. In 
presence of Rashba interaction ground states have a dif- 
ferent structure and their polarization varies with the 
density. However the ground state energy is always lower 
than in absence of spin-orbit interaction over the whole 
density range. 

In Fig. [S]the pair correlation function g(r) decomposed 
into triplet and singlet components is reported. Triplet 
and singlet components were obtained collecting the fol- 
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FIG. 3: (color online) Energy values at rs = 10 A = 0.05 as in 
the figure above. DMC results are represented by triangles, 
HF by squares and results in absence of Coulomb interaction 
by circles. 



The results show a tendency for the triplet and singlet 
first peaks to become smaller and closer to eachother 
when spin-orbit interaction is switched on. This sug- 
gests a decrease in the antiferromagnetic character of the 
electrons, probably due to the induced spin polarization. 
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FIG. 4: Difference between energies without and in presence 
of Rashba interaction (A = 1). Data were taken from^^. 



lowing quantities: 



g,(.) = jv^<^^''^f;;'^,J'^^^ (30) 



{ipTlSjrij - r)(Ti ■ (Tj\ipT) 

{ipT\i>T) 



(31) 



where A'' is a normalization factor, and composing them 
according to the triplet and singlet spin projectors'^: 



(32) 
(33) 



FIG. 5: (color online) Pair correlation function g(r) decom- 
posed into singlet (circles) and triplet (squares) channels for 
T's = 5 with (closed symbols) and without (open symbols) 
spin-orbit interaction. For better comparison the curves are 
all normalized to give g{r) — » 1 for r — > oo 



IV. CONCLUSIONS 

We have implemented a DMC algorithm for treating 
Rashba spin-orbit interaction in the electron gas. We 
have calculated the equation of state of 2D electron gas 
in the presence of Rashba interaction. We computed the 
energy per particle as a function of the strength A of 
the Rashba potential, and the polarization ^ showing 
system polarization and comparing the results with an- 
alytical HF results. Our work not only gives results for 
the 2D electron gas, but also provides a good test for the 
algorithm, opening new possibilities for implementation 
on other systems, in particular quantum dots or quan- 
tum wires involving both the Rashba and the Dresselhaus 
potentials!. A similar treatment can also be extended to 
finite systems such as atoms and molecules. 
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APPENDIX A: TABLES WITH RESULTS 

In the following tables the Diffusion Monte Carlo re- 
sults are reported together with corresponding Hartree- 
Fock energies and results for the system without the 
Coulomb interaction with the Rashba interaction only. 
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TABLE L Energies of the system at = 1 and A = 0.1, 0.2, 
and 0.5 for various polarizations ^ as defined in Eq. (|28p . The 
DMC result is computed with Coulomb and Rashba interac- 
tion, Enashba IS the energy given by the Rashba Hamiltonian 
(|25p . and in the last column the Hartree-Fock result is re- 
ported. All the energies are expressed in Ry. 



A = 0.1 




TV- 


N+ 


EoMC 


Enaahba 


Ehf 


0.690 


9 


49 


0.243(4) 


1.67 


0.244 


0.414 


17 


41 


-0.128(3) 


1.29 


0.009 


0.276 


21 


37 


-0.259(3) 


1.15 


-0.081 


0.138 


25 


33 


-0.354(2) 


1.06 


-0.151 


0.069 


27 


31 


-0.388(4) 


1.024 


-0.178 


0.0 


29 


29 


-0.414(3) 


1.00 


-0.200 


-0.069 


31 


27 


-0.427(4) 


0.985 


-0.217 


-0.138 


33 


25 


-0.432(3) 


0.98 


-0.229 


-0.207 


35 


23 


-0.428(3) 


0.984 


-0.235 


-0.276 


37 


21 


-0.416(5) 


0.998 


-0.237 


-0.414 


41 


17 


-0.360(4) 


1.06 


-0.223 


-0.690 


49 


9 


-0.136(4) 


1.29 


-0.137 


A = 0.2 




TV- 


N+ 


Edmc 


E flashba 


Ehf 


0.690 


9 


49 


0.437(4) 


1.86 


0.435 


0.414 


17 


41 


-0.008(4) 


1.40 


0.125 


0.276 


21 


37 


-0.184(5) 


1.23 


-0.022 


0.138 


25 


33 


-0.306(5) 


1.1 


-0.112 


0.0 


29 


29 


-0.411(4) 


1.00 


-0.200 


-0.138 


33 


25 


-0.469(3) 


0.94 


-0.268 


-0.276 


37 


21 


-0.514(5) 


0.92 


-0.314 


-0.414 


41 


17 


-0.474(5) 


0.94 


-0.340 


-0.690 


49 


9 


-0.322(4) 


1.09 


-0.328 


A = 0.5 




iV- 


N+ 


Edmc 


Ejlashba 


Ehf 


0.690 


9 


49 


1.013(7) 


2.43 


1.008 


0.276 


21 


37 


0.053(6) 


1.47 


0.230 


0.0 


29 


29 


-0.410(6) 


1.00 


-0.200 


-0.207 


35 


23 


-0.660(3) 


0.75 


-0.469 


-0.414 


41 


17 


-0.821(5) 


0.59 


-0.688 


-0.690 


49 


9 


-0.888(8) 


0.52 


-0.900 


-1.0 


58 





-0.741(7) 


0.67 


-1.031 



9 



TABLE II: Energies of the system at = 5 and A = 0.02, 
and 0.1 for various polarizations ^ as defined in Eq. (|28p . The 
DMC result is computed with Coulomb and Rashba interac- 
tion, Enashba IS the energy given by the Rashba Hamiltonian 
(|25p . and in the last column the Hartree-Fock result is re- 
ported. All the energies are expressed in Ry. 



A = 0.02 




iV_ 


N+ 


Edmc 




Ehf 


0.690 


9 


49 


-0.2719(5) 


0.0667 


-0.2178 


0.414 


17 


41 


-0.2861(5) 


0.0515 


-0.2042 


0.138 


25 


33 


-0.2947(5) 


0.0423 


-0.1995 


0.0 


29 


29 


-0.2971(5) 


0.0400 


-0.2001 


-0.138 


33 


25 


-0.2978(4) 


0.0392 


-0.2026 


-0.414 


41 


17 


-0.2951(5) 


0.0422 


-0.2125 


-0.690 


49 


9 


-0.2867(4) 


0.0514 


-0.2330 


-1.0 


58 





-0.2748(5) 


0.0693 


-0.2702 


A = 0.1 




iV_ 


N+ 


Ebmc 


E Rashba 


Ehf 


0.690 


9 


49 


-0.2403(5) 


0.0972 


-0.1872 


0.414 


17 


41 


-0.2668(4) 


0.0701 


-0.1856 


0.138 


25 


33 


-0.2879(4) 


0.0486 


-0.1932 


0.0 


29 


29 


-0.2964(4) 


0.0400 


-0.2001 


-0.138 


33 


25 


-0.3033(3) 


0.0330 


-0.2088 


-0.414 


41 


17 


-0.3127(3) 


0.0236 


-0.2321 


-0.690 


49 


9 


-0.3155(5) 


0.0209 


-0.2635 


-1.0 


58 





-0.3099(5) 


0.0267 


-0.3129 



10 



TABLE III: Energies of the system at = 10 and A = 0.02, 
0.05 and 0.1 for various polarizations ^ as defined in Eq. (|28|l . 
The DMC result is computed with Coulomb and Rashba in- 
teraction, Enashba is the energy given by the Rashba Hamil- 
tonian (|25p . and in the last column the Hartree-Fock result 
is reported. All the energies are expressed in Ry. 



A = 0.02 




iV_ 


N+ 


Edmc 


E flashba 


Ehf 


0.690 


9 


49 


-0.1621(5) 


0.0186 


-0.1236 


0.414 


17 


41 


-0.1660(2) 


0.0140 


-0.1138 


0.138 


25 


33 


-0.1688(3) 


0.0110 


-0.1099 


0.0 


29 


29 


-0.1694(2) 


0.0100 


-0.1100 


-0.138 


33 


25 


-0.1698(2) 


0.0094 


-0.1149 


-0.414 


41 


17 


-0.1700(2) 


0.0094 


-0.1184 


-0.690 


49 


9 


-0.1699(2) 


0.0109 


-0.1313 


-1.0 


58 





-0.1684(3) 


0.0147 


-0.1551 


A = 0.05 




iV_ 


N+ 


EoMC 


E flashba 


Ehf 


0.690 


9 


49 


-0.1561(4) 


0.0243 


-0.1179 


0.414 


17 


41 


-0.1625(3) 


0.0175 


-0.1103 


0.138 


25 


33 


-0.1669(3) 


0.0122 


-0.1088 


0.0 


29 


29 


-0.1690(2) 


0.0100 


-0.1100 


-0.138 


33 


25 


-0.1714(2) 


0.0082 


-0.1127 


-0.414 


41 


17 


-0.1736(2) 


0.0059 


-0.1219 


-0.690 


49 


9 


-0.1746(2) 


0.0052 


-0.1370 


-1.0 


58 





-0.1742(2) 


0.0067 


-0.1631 


A = 0.1 




iV_ 


N+ 


Edmc 


Ejlashba 


Ehf 


0.690 


9 


49 


-0.1459(4) 


0.0338 


-0.1084 


0.414 


17 


41 


-0.1565(3) 


0.0233 


-0.1045 


0.138 


25 


33 


-0.1656(3) 


0.0141 


-0.1068 


0.0 


29 


29 


-0.1696(2) 


0.0100 


-0.1100 


-0.138 


33 


25 


-0.1733(2) 


0.0063 


-0.1146 


-0.414 


41 


17 


-0.1792(2) 


0.0001 


-0.1277 


-0.690 


49 


9 


-0.1837(2) 


-0.0043 


-0.1465 


-1.0 


58 





-0.1861(2) 


-0.0067 


-0.1764 



11 



TABLE IV: Energies of the system at = 20 and A = 0.01, 
0.02 and 0.1 for various polarizations ^ as defined in Eq. (|28|l . 
The DMC result is computed with Coulomb and Rashba in- 
teraction, Enashba is the energy given by the Rashba Hamil- 
tonian (|25p . and in the last column the Hartree-Fock result 
is reported. All the energies are expressed in Ry. 



A = 0.01 




iV_ 


N+ 


Ebmc 


E flashba 


Ehf 


0.690 


9 


49 


-0.0901(3) 


0.0046 


-0.0665 


0.414 


17 


41 


-0.0912(2) 


0.0035 


-0.0604 


0.138 


25 


33 


-0.0919(2) 


0.0027 


-0.0577 


0.0 


29 


29 


-0.0921(2) 


0.0025 


-0.0575 


-0.138 


33 


25 


-0.0923(2) 


0.0024 


-0.0581 


-0.414 


41 


17 


-0.0923(3) 


0.0023 


-0.0616 


-0.690 


49 


9 


-0.0920(3) 


0.0017 


-0.0684 


-1.0 


58 





-0.0912(4) 


0.0037 


-0.0812 


A = 0.02 




iV_ 


N+ 


EoMC 


Ejlashba 


Ehf 


0.690 


9 


49 


-0.0892(2) 


0.0056 


-0.0655 


0.414 


17 


41 


-0.0906(2) 


0.0041 


-0.0598 


0.138 


25 


33 


-0.0918(3) 


0.0029 


-0.0575 


0.0 


29 


29 


-0.0922(2) 


0.0025 


-0.0575 


-0.138 


33 


25 


-0.0924(2) 


0.0022 


-0.0583 


-0.414 


41 


17 


-0.0928(2) 


0.0018 


-0.0622 


-0.690 


49 


9 


-0.0928(3) 


0.0018 


-0.0693 


-1.0 


58 





-0.0924(5) 


0.0023 


-0.0825 


A = 0.1 




iV_ 


N+ 


Edmc 


ERashba 


Ehf 


0.690 


9 


49 


-0.0811(3) 


0.0132 


-0.0579 


0.414 


17 


41 


-0.0856(2) 


0.0087 


-0.0552 


0.0 


29 


29 


-0.0918(2) 


0.0045 


-0.0560 


-0.138 


33 


25 


-0.0937(2) 


0.0006 


-0.0599 


-0.414 


41 


17 


-0.0972(2) 


-0.0029 


-0.0668 


-0.690 


49 


9 


-0.1001(2) 


-0.0058 


-0.0770 


-1.0 


58 





-0.1026(2) 


-0.0083 


-0.0932 



